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Abstract. We construct analytic solutions to the Euler equations with an interface between two 
fluids, extending work of Duchon and Robert. 

1. Introduction 

We consider the interface problem for two perfect incompressible fluids in R 2 under the influence of 
gravity. Each of the two fluids occupies one of the two connected regions f2±(t) in the complement 
of a moving interface parametrized by a curve y(t, x\): 

£(i) = {(xi,x 2 ) : x 2 = y{t, xi)}. 

The relevant equations of motion are given by the two-dimensional Euler system: 



(1) 



Vt + V • {p±v <£> u) + Vp = p±g in fi± 
div v = in 



where v = (vi, v 2 ) is the velocity of the fluid, p+ =/= p- are the densities of the two fluids, and g > 

is the gravitational acceleration. We also introduce the Atwood number a = p+ ~ p ~ ■ 
to p++p- 

The assumption of incomprcssibility implies that the vorticity fi := V x u is a distribution supported 
on S. We are thus interested in studying vortex sheet solutions to the system (JXJ) , i.e. solutions 
where Q is a measure supported on E. In |DR| . the authors construct global solutions to the vortex 
sheet problem in a homogeneous fluid, i.e. in the case p+ = p_, with small initial data. The object 
of the present paper is to extend their result to the the case a =/= 0. 

2. Equations for vortex sheets 

We begin by presenting the equations for vortex sheets. These arc formulated in terms of the 
vorticity density ui, defined by: 

(SI,® = [ Cj^x^y^x))^ 
Jm. 

for all <fi G Cj?°. The components of the velocity field can be expressed in terms of uj by the 
Birkhoff-Rott integrals: 

, +s 1 f yfat) - y(x',t) , 

^ = -^1 (x-x>r + (y(x ; t)-y(x>,t))^ X > t)dX 

V2M = i P ' V - / (x - x> ? + (y^t) - y(x>, t)f t] dX '- 

l 



(2) 
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Sulem, Sulem, Bardos and Frisch |SSBF| have derived the equations of the motion of a vortex 
sheet: 

d t y + viy x = v 2 



(3) 



dtiljtb - a(vi + v 2 y x )) + d x {vi(^u> - a(yi + v 2 y x )) 

u 2 \v\ 2 



Following [PR] , we begin by expressing the system ([2]), Q in a form that isolates its linear part. 
Introducing the function u), defined by 

ui = 2(l + a;), 

and differentiating the first equation with respect to x, we can rewrite (|2|). (|3|) as 

d t y x -Kw = N l 

(4) < d t u + aHdtVx - Ay x - ad x u) + agy x = N 2 

Vx(t = 0)=y x o 

where N± = N\{y x ,uj) and ^2 = N 2 (y x ,uj) arc nonlinear terms whose exact form will be given in 
section 4 (see (|12[) below). A is the pseudo-differential operator with symbol |£|, and H is 

the Hilbert transform, with symbol — The equations Q form the vortex sheet system with 

which we shall be concerned. 

Before stating our result, let us recall the definition of the spaces introduced in |DRj . which will be 
used to carry out the argument. 

Definition 2.1. We denote by Bo the space of function whose Fourier transforms are bounded 
measures, with the norm 

r d\u\. 



For p > 0, the space B p C Bq is defined by its norm: 



d|e'KlTi(0|. 



For a > 0, B a is the subspace of C® (R; Bq) defined by the norm 

ll^llsc = inf { s : M a positive bounded measure, \e at ^u(t)\ < fi for all t > 0}. 
In addition, \n\ a denotes the measure achieving the infimum. 

The Paley- Wiener theorem ensures that functions in B p extend analytically to the strip {\?Sz\ < p] 
of the complex plane. In addition, B p is a convolution algebra. 

The main result of the present work is the following: 

Theorem 1. Let < \a\ < 1 be fixed with ag < 0. For some e\ > 0, for any yo x G Bq with norm 
< ex, there exists a > and a solution (y x ,uj) in B a x B a to the system Q). Such a solution is 
unique in some ball of Bq x Bq, and y x (t), uj(t) converge to zero uniformly as t goes to infinity. 
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Given ag > 0, and any fixed time T > then there are a > 0, e 2 {T) > such that for any yo x € Bq 
with \\yx\\Bo ^ e 2; there exists a pair (y x ,u>) in B a x B a which solves the system for < t < T . 

We begin the proof of the theorem in the case ag < in the next section, where we further reduce 
the system ((4|) to a form suitable to the application of a contraction mapping argument in B a . In 
section [5l we explain how to adapt the estimates in |DR| to the nonlinear terms we obtain in the 
case a ^ 0. Finally, in section 6 we show how to address the case ag > 0, and establish the second 
part of the theorem. 



3. The linear system 

It will be convenient to rewrite (j4]) in terms of suitable combinations of y x and u>. We consider the 
case Ni = N2 = 0. Substituting the first equation into the second, we find, assuming decay for 
t — > 00: 

(5) Q = (ia ill ± v 7 ! - a 2 - offllp 1 )^ 
or, in terms of the differential operator A, 

(6) uj = {-aH ± yjl-a 2 -agk- l )y x . 
This suggests rewriting the system (j4]) as 



(7) 



d t (± VI - a 2 - agA-iy? ) = ± y/l - a 2 - agk^ku ± y/l - a 2 - agA- 1 ^ 
d t (u + aHy x ) = Ay x + ad x Lo - agy x + N 2 . 



Assuming Ni = N2 = and adding the two equations above, we obtain the following pair of 

equations: 

(8) 

d t {u + aHy x ± ^\ - a 2 - agA- 1 y x ) = (ad x ±\/l - a 2 - agA- l A){io+aHy x ±^Jl - a 2 - agA- 1 y x ). 
To study the linear equations ([7]), we introduce the notations: 



A Q . g = y/1 - a 2 - agA~\ 
X±=uj + aHy x ± A a , g y x , 
S± = e *( Q ^±A a , g A) ; 

I+h{t) = / S+(t-s)h(s)ds, 



and the operators: 



I~h{t) = / S-(t-s)h(s)ds, 
J a 

I h = l + h(0). 

We can formally write solutions to ([7]) as 

X+(t) = S+(t)X+(0) - I+(N 2 + A a _ g Ni)(t) + S+(t)I Q (N 2 + A , fl JVi), 
X-(t) = 5_(<)A_(0) + r(N 2 - A^gNtXt). 
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Setting .X+(0) = — io(^2 + A ajg 7Vi) puts X+(t) in Bo, and we have the relation 
(9) w(0) = -aHy x0 - A a , g y x0 - I (N 2 + A^JVi). 

Taking equation ([9]) into account, we express y x and uj in terms of X + and X_. We obtain the 
following system of equations: 



(10) 



y x = S-(t)y x0 - 2^(- 5_(t)J (JV a + A Q)9 iVx) 
+ /+(iV 2 + A„, fl JVi) + r(N 2 - A^JVi)), 
w + affife = S-(t)(-A aifl ifco) - -(S-(t)I (N 2 + A 0iS iVi) 



- I+(JV 2 + A„, fl JVi) - J-(JV 2 - Ao^JVi)), 
k 2/ x (i = 0) = 2/^o- 
We will apply a fixed point argument to this last system. 

Before analyzing the nonlinear terms, let us derive two simple estimates for the operators J+ and 
7_. 

Lemma 3.1. Suppose ag < 0. There exists ao > such that for < a < ao o,nd any F £ B a , we 
have the estimates: 

l^i^U < C(a)\F\ a , 
[I^Ftlc < C{a)\F\ a . 

Proof. We will prove the estimates for the operator I~ , the estimates for I + follow similarly. 









< 


/' 






Jo 




< 


\F\ a / 






■-/(i 




\F\ a — 
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Aa, 9 (0l 



(1-e 



t|e|(o-Aa,„)l 



When ag < 0, the last line is uniformly bounded in t if a is sufficiently small. 
To derive the analogous bound for I~ F t , we integrate by parts in s: 



e at ^\I-F t \ 



e «l€l* / e (t- S )(mS-A a , s (S)|«|) 5 



< 



s=t 



s=0 

3 (t_ 8 )( io e-A(i)i«i) • 



< e a l e l*|^,t)| + e ( a - A «-« tt Wl c l*|F(f,0)| 
+ I^U f |£|(|o| + A Otff (0)e( t ->(° t - A -«Wl c l d S 



< C(a)\F\ a . 
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4. The Structure of the Non-Linear Terms 

We consider the structure of the non- linear terms Ni(y x ,uj) and N 2 (y x ,uj) appearing in the system 
(H|). From the equations ([3]), (jl]) we find: 

(11) N 1 (y x ,oj) = (-Viy x + (v2-Hoj)) x 

and 

N 2 {y x ,u) = ((avx + aHy x ) + av 2 y x ) t 



(12) + ( UJVi + (v\ + Hy x ) — av 1 — av\Viy x — a 



2(1 + 1,1) "V 2 2 



The subscripts x and t denote partial differentiation in the variables x and t, respectively. We have 
grouped together terms with their linear part. 

We note that Ni(y x ,uj) = F x has no dependence on the parameter a. N2{y x ,w) depends on a and 
is different from the function G(y x ,uj) x appearing in jDRj . p. 217. However, we can still write 
N2{y x , w) in the form 

N 2 (y x ,u) = Gi(y x ,ui)t + G 2 (y x ,w) x . 
Suppressing the dependence on t, we define 



p = p(x,x') 

In terms of p, the equations ((2]) then become: 



y(x) - y(x') 



(13) 



1 f p 1 + u^Q , 

fiW = P- v - / -i , 2 — , 

7r J 1 + p a; — a; 

v 2 (x) = - P-v. / 2 — dx . 

1" J 1 + P X — X' 



We recall the following singular integral representations for the operators H and A: 

» / \ 1 f u(x) - u(x') , 1 f u(x') , 

Au(x) = -p.v. / — - 1 7To~~^ dx and = -p.v. / , dx . 

7T J (X — X Y 7T J X — X 

We will formally expand F(y x ,uj), G\{y x ,ui) and G2{y x ,oj) in terms of the sequence of singular 
integral operators {Tj}j>i defined by 

rp ( \ ( \ 1 [ ( y(x)-y{x') \ j u{x') A , 

Tj(y x )u{x) := - p.v. / - dx . 

7T J \ X — X'/X — X' 

For F(y x ,uj), we have, just as in |DRj : 

1 /" / 1 \l+w(x') , . 1 /• p l+w(x') , . 

7T J X — X 7T J 1 + X — X 



and thus, 



F(y x ,ui) = ^e j T j {y x )(l + u)+y x ^e' j T j (y x )(l+ui), 

3=2 3=1 
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where €j, e'j equals — 1, or 1. Note that this expression is at least quadratic in y x and u>. 

We now want to show that G\{y x ,ui) and G 2 {y x ,bS) have a similar decomposition. From equation 
(|12[) . the evolution equation for u, we see that one term in G\{y x ,Lo) that we need to consider is 
av\ + aHy x . Using the expression for v\(x) in (|13j). we have 

oo 

(14) avi +aHy x = -aT^y^u ■ </ ^ < // ,•://. :■; 1 

3=3 

Similarly, using the expression for 1*2(2;), the other term in G\{y x ,ui) satisfies 

^ 00 
av 2 y x = -y x Hu; + y x ^ e j T j (y x )(l + uj). 

7T * — ' 

We see that, just as for F(y x ,u>), every term in the expression for G\ (y x ,u>) either contains a factor 
of Tj(y x )(l + uj) for some j > 2, or else at least two factors of y x and ui, i.e. G\ is quadratic in y x 
and to around (0,0). 

Turning to G2 (y x ,uj), each term can either be written in the form above, or else is of the following 
type: 

/ OO v f OO 

(15) aviv 2 y x = ay x f ^ e j T j (y x )(l + uj) J f Huj + ^ <':r,(j/, n 1 + 10) 

Each term in G2{y x ,oS) either contains a factor of Tj(y x ){l + w) for some j > 2 or else contains a 
factor of T\(y x ) together with at least one other factor of y^, w or T\{y x ). 

5. The nonlinear estimate 
The key estimate of Duchon and Robert's work is 

Lemma 5.1 ([DR], Section 4). Lety\ x , y 2x , U)i, w 2 be elements ofB p . Assume \yi X \ p < \i, \^%\ P < A* 
(i = 1, 2) and \yi x — y 2x \ p < v, \u)\ — u 2 \ p < v with \i and v positive bounded measures, J d[i < 1. 
Then, 

\F{yi x ,uji) - F{y 2x ,uj 2 )\ p < A(fi) * v, 
where A is a continuous map from the open unit ball of A4 + (the set of bounded positive measures) 
into M + , with A(0) = 0. 

We show in this section that the analogous inequality for Gi(y x ,cj) and G 2 (y x ,u>) holds. To prove 
the above lemma, Duchon and Robert expressed F{y\ x ,uii) — F(y 2x ,uj 2 ) in terms of the operators 
Rk(vix, ■ ■ ■ , Vkx), defined by 

1 f 1 

Rk(yix, ■ ■ -,ykx)tt(x) := - p.v. / -(pi ■ ■ ■ p k ) x Q,(x ) dx , 
where pi(x,x') = (yi(x) — yi{x'))/{x — x'). By abuse of notation, we also set 

Rk{y x ) ■■= Rkijjx, ■ ■ -,y x )- 

They then used the following proposition: 
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Proposition 5.1. Let y% x , . . . , y^x be elements of B p . Then Rk{y\ x , ■ ■ ■ , Vkx) maps Bq into B p , 
and \Rk(yix, ■ ■ -,ykx)^\p < *•••* \ykx\p * l^l-p, where |£l|_ p stands for exp(-7rp|£|)|f2| . 

In particular, this allows us to deduce that 

(16) \Tj(v*M P < (i + 2j)l»«i; J *|fiU 

and 

(17) iT^y^t - T 3 {y 2x )n 2 \ p < c(j)(ti<i-V + ^) * „, 
where c(j) := 2j 2 + 3j + 1 and fj,*^ := \x * ■ ■ ■ * fx (j times). 

For any terms in G\(y x ,uj) and G 2 (y x ,uj), of the form found in equation (|14p. we can proceed in 
an identical fashion to Duchon- Robert to establish Lemma I5TT1 



We are left with terms of the form in equation (fT5j) . We split this term into two parts. Firstly, 
consider 



*l(j/ x ,w) := ay x ■ Hu ■ ( ^ ejTj(y x )(l +w) V 
V j= i / 



Since H is an isometry in B p , we can use equations (|16[) and f|l T[) to show that 

|*l(yix,Wl) - *l(j/2a:,W2)|p < ^l(^) * ^ 

for a continuous function Ai satisfying A\(0) = 0. 
We also need to consider 



oo 



^2(y x ,uj) ■= ay x ■ ( ^ejTj(y x )(l +w) j ■ ( ^ e^Tj(y x )(l + w) 

Using equations f| 16[) and (|17p again, we have the same inequality for ^ for some continuous 
function A 2 satisfying A 2 (0) = 0. 

To summarize, for \y lx \ p , \u>i\ p < fi and \y\ x - y^x\ P , \wi - uj 2 \ p < v, we have 

(18) ^(yix.wi) - Fdfts,^)^ < * 

(19) \G 1 (y lx ,u>i) + G 2 (yi x ,u>i) - Gi(y 2xi uj 2 ) - G 2 {y 2xi uj 2 )\ p < A(p) * u, 

for a continuous function A satisfying A(0) = 0. 

Using ([15]). (fl9|) . one can now easily construct solutions to (|T0|) with prescribed initial data y x G -Bo 
and w(0) given by (|9]) by a contraction argument in B a x B a around the linear solution 

Ux,lm = S-(t)y xo 
Win = -aHS-(t)y Xo - S_(i)(A O)g y x0 ). 
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6. The case ag > 

In this section, we prove the second part of Theorem [TJ When ag > 0, we cannot reproduce 
the estimates in Lemma 13.11 because the linear semigroup does not dampen small frequencies in 
time. 

We now express the solutions of © in a form adapted to the linear analysis we wish to carry out. 
Let x be a smooth cut-off function such that x = 1 f° r 1*1 < 1 an d X = for t > 2. Denote 
Xt(0 = X(-/T). 



We begin by setting 



I + h(t) = / S+(t-s) X T(s)h(s)ds, 



rh(t)= S-(t-s)xT(s)h(s)ds 
J a 

I h = (l + )h(0). 



For t < T, we have, 



d t I + h(t) = h(t) + (ad x + ^Jl - a 2 - agK- 1 K)I + h{t), 
and the corresponding equation holds for I~ . Thus, if we let 

X_(t) = S^(t)X_{0)+r(N 2 - A at gNi)(t), 

and 

X+(t) = S+(t)X+(0) - I+(N 2 + A^mXt) + S+(t)I (N 2 + A^iVx), 
X + and X_ satisfy equations ([7]) for i < T. 

As before, we have only prescribed initial data for y x and not for u>. This means that we still have 
one degree of freedom, and we use this to ensure that X±(t) are in Bo- Namely, we set 

X+(0) = -J (JV a + A a , 9 7V 1 ). 

Thus, 

X_(0) = -2A OiS y x0 - J (JV 2 + A„ )4 ,JVi), 

and 

w(0) = -aii^o - A OiS y x0 - / (^ 2 + A 0iS iV~i). 

Since the difference of any two solutions will satisfy the homogenous equations, the solutions above 
are unique in Bo x Bo up to adding a constant to w(0). Choosing w(0) to have mean zero selects a 
unique solution. 

Rearranging the above equations, we find that 

2y x = 2S_(i)y x0 + A~* (S_(*)I (JV a ) - ^(JV a )(t) 

- 7+(iV 2 )(i)) + 5T_(t)/ (JVi) + I~(Ni)(t) - I+iNJit), 
20 < _ 

2(w + offy x ) = -2S_(i)A o , s2 / x0 - S-(t)I (N 3 + A a , fl JVi) + /-(JV 2 - A„, fl JV 1 )(t) 

-7+(iV 2 + A 0)S iVi)(i). 



DUCHON-ROBERT SOLUTIONS FOR A TWO-FLUID INTERFACE 



9 



A solution (y x ,u>) <G B a x B a of these equations is a local solution of the vortex sheet equations for 
t < T. 

If we choose yo x <G Bo, then S-(t)y x o € B a for small a > 0. We can thus hope to apply a fixed 
point analysis in B a x B a around (S-(t)y X Q, ~S-(t)A a ^ g y x o) to the above formulation. 

The substitute for Lemma T3. II is the following: 

Lemma 6.1. Write, as in section 4- 

N 1= F X 

N 2 = (G l ) t + (G 2 ) x . 

There exists ao > such that for < a < a^, the nonlinear terms on the right side of A20\) are 
bounded by 

e CT C{a)(\F\ a + \Gi\ a + \G 2 \ a ). 

Proof. Let us first consider the high frequencies 

|£| > 3o fl /(l - a 2 ). 

For this range of |£|, A a ff is real and both bounded and bounded away from zero. We can bound 
the above operators as in Lemma 13. II 

We also consider the low frequencies, |£| < 3a<7/(l — a 2 ). There are various terms to consider: 
(1) f-(F x )(t) and I + (F x )(t): For I+(F x )(t), we have: 

/DO 

<TC(a)\F\ a . 



We have also used |£| is bounded in the region of interest. A similar procedure leads to 

\I-F x \ a <e CT C(a)\F\ a . 

(2) I~ ((Gi)t)(t) and I + ((Gi)t)(t): since ia£± A aiS |^| is bounded uniformly for |£| bounded, we 
can integrate by parts and bound the integrand uniformly to conclude the same inequalities 
as in (1). 

(3) A a ^ g I~ (F x )(t) and A a , g I + (F x )(t): since £A aj9 remains bounded for small £, we can bound 
the integral uniformly to conclude the same inequalities. 



(4) A-i (M0V(G 2 )*) - 7-((G 2 )*)(*) - 7+((G 2 ) x )(t)J. This case requires more work. We 
can express this term as 

(21) A-* (J (S-(jk)S+(-s) -S-(t- s)) X T(s)(G 2 ) x (s)ds^ 

+ A-] ( [°° (S-(t)S+(-s) - S+(t - s)) X T(s)(G 2 ) x (s)ds) 
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The multiplier A a g vanishes at |£| = yz^r- Away from this frequency, we can proceed as 
in the previous cases to obtain the desired bounds. Near |£| = j^z, we have the following 
bound: 



\F((S-(t)S + (-s)-S_(t-s))f) (0| 

< e t(ta£-A , a |£|) e s(fo£-A ,g|£|) _ e (*-s)(ia£-A a , 9 |£|) 

<qA a , a |e sA -l^|/(e)|. 



1/(01 



Inserting this expression into the first term in (|2ip. we obtain the required bound. The 
second term in (|2T|) is bounded analogously. 



(5) A"! (V(i)io((Gi)t) - J-((Gi) t )(t) - J+((Gi) t )(*)V This term is bounded by the same 
procedure as the term in the previous item. 

□ 

We see that for T ~ logi, the bounds in Lemma 16. 1 1 are of order e _1 . Using these estimates, we 
can repeat the fixed point argument as in the case ag < and obtain a local solution in a ball in 
B a x B a for f < T. 

We can obtain a somewhat improved time of existence by introducing the following norm. 

Definition 6.1. For a > 0, £> Q is £/ie subspace of Ct(R; .Bo) defined by the norm 

\\u\\g = inf{||u||s : /i a positive bounded measure, f(£,t)\u(t)\ < /i for all t > 0}, 

where /(£, i) = e Q *'^ /or £ < 1 and f(£,i) = e Q '^ otherwise. The measure achieving the infimum is 
denoted \u\a ■ 



Working in the space B a , we can improve the bounds in Lemma 16.11 to estimates of order 

TC(a)(\F\ Sa + \ Gl \ Sa + \G 2 \ Sa ). 

To use the fixed point argument, we need to check that the analog of Proposition 15.11 holds for the 
space B a . We have the following: 

Proposition 6.1. Let y\ x , . . . , y^ x be elements of B p . Then Rk(yix, ■ ■ ■ , ykx) enjoys the bound 

\Rk{y\x,---,ykx)&\* < C\y\x\& * • •• * \vkx\z * Mb „■ 



The proof follows readily from Proposition ^. II Thus, the improved dependence on T of the bounds 
above allows us to extend the time of existence to T ~ r 1 for initial data of size e in B . 
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